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The Improvement Program in Nonrelativistic Lattice QCD 

Colin J. Morningstar a * 

a Stanford Linear Accelerator Center, Stanford University, Stanford, California 94309 

Progress in the improvement program in nonrelativistic lattice QCD is outlined. The leading radiative cor- 
rections to the heavy-quark mass renormalization, energy shift, and two important kinetic coupling coefficients 
are presented. The reliability of tadpole-improved perturbation theory in determining the energy shift and mass 
renormalization is demonstrated. 
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1. INTRODUCTION 

The action in nonrelativistic lattice QCD [1,2] 
(NRQCD) must be engineered to reproduce the 
action of continuum QCD at low energies. The 
starting point in this process is a discretized ver- 
sion of the Schrodinger action. Interactions must 
then be added to systematically correct for rel- 
ativity and finite lattice spacing errors. Lastly, 
the coupling strengths of these added interac- 
tions must be determined. One possibility is to 
treat these new couplings as adjustable param- 
eters and tune them to fit certain experimental 
data; however, this tuning significantly reduces 
the predictive power of NRQCD simulations and 
is very costly and difficult. A better alternative 
is to compute the new couplings in terms of the 
fundamental QCD coupling a s and the heavy- 
quark mass M using perturbation theory. This 
is done by evaluating various scattering ampli- 
tudes in both QCD and lattice NRQCD and ad- 
justing the couplings until these amplitudes agree 
at low energies. Since the role of these couplings 
is to compensate for neglected low-energy effects 
from highly-ultraviolet QCD processes, one ex- 
pects that they may be computed to a good ap- 
proximation using perturbation theory, provided 
that M is large enough and the lattice spacing 
a is small enough. This entire engineering pro- 
cedure constitutes the improvement program in 
NRQCD. 

In this talk, I report on progress being made 
in this improvement program. I present results 



for the leading radiative corrections to the mass 
renormalization, energy shift, and two important 
kinetic coupling coefficients. This work is an ex- 
tension of previous calculations [3,4] and is being 
carried out in conjunction with ongoing simula- 
tions as part of the NRQCD collaboration [5-7]. 

2. LATTICE NRQCD 

The lattice NRQCD action used here is written 
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where n is a positive integer, tp(x) is the heavy- 
quark field, Ufj,{x) is the link variable representing 
the gauge field along the link between sites x and 
x + ae^, and uo is a mean-field parameter [8] de- 
fined in terms of the mean plaquette. The remain- 
ing operators are given by Ho = — A^ 2 ) /(2M) and 
5H = Ej=i Cj Vj , where the Vj denote various 
interactions constructed from covariant, tadpole- 
improved difference operators and the cloverleaf 
chromoelectric and magnetic fields. Two of these 
operators are given by 
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where a 2 "A( 2m ) = ELiK \U k + Ul]-2] m . The 
parameter n stabilizes the evolution of the quark 




Figure 1. The energy shift parameter A, the heavy-quark mass renormalization parameter B, and the 

(2) (2) 

kinetic coupling coefficients c\ and c 2 , against the bare heavy-quark mass aM. For aM > 3, the 
stability parameter n is set to unity; for 1.5 < aM < 3, n = 2 is used; for 1 < aM < 1.5, n = 3 is used; 
and for 0.5 < aM < 1, n = 6 is used. 



Green's function when n > 3/(Ma), approxi- 
mately. The coupling coefficients Cj are functions 
of a s and aM; at tree level, their values are unity. 

3. HEAVY-QUARK SELF-ENERGY 

The 0(a s ) contribution to the heavy-quark 
self-energy has been determined. The small-p 
representation of the self-energy may be written: 

2 2 / 2\2 4 

aE(p) « oJ/oHi/iH — +/ 2 H^A 
\fj s^joy j^jiy ) 2 Ma J V ' 8M 2 a 2 

+ /sW (X>> 4 ) + •••}■ 

i=i 

where = —ip^a and f m {w) = J^fco ^™ The 
on-mass-shell quark satisfies a dispersion relation 

( ^ p2 ( p2 ) 2 * r u 



where M r = Z m M is the renormalized mass with 
Z m = 1 + a s B, B = Q^ + Q^ + Q^, and 

aMp) = ^ 0) +4 0) (Ef> 4 ) + ^^2^. 

i=i 

+ 2n( 2 >+n<°>('2+ ] i-)+2n('>+n< >. 

The sole effect of the 0(a s ) corrections in QCD 
is to renormalize the quark field and mass. If lat- 
tice NRQCD is to reproduce the low-energy phys- 
ical predictions of full QCD, then the 0(a s ) cor- 
rections to the heavy-quark propagator in lattice 
NRQCD must also do no more than renormal- 
ize the heavy-quark field and mass at low mo- 
menta. This will be true if 5w(p) = 0; that is, if 
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Table 1 



Comparison of the T mass expected from perturbation theory using Eq. 2 and the mass M. 

obtained by fitting simulation results to the dispersion relation in Eq. 1. Values for ay(q*) are obtained 
by fitting simulation measurements of the mean plaquette to its perturbative expansion [5]. 
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constant term Q 



(o) 



= 0. Alternatively, since the 
represents an overall energy 



shift, one could require only Q^ 0) = Q w = and 
then simply shift the energies obtained in sim- 
ulations by an amount —a s A/a for each heavy 



quark, where A = — Q 



(o) 



Setting fi^ 







then determines the kinetic coupling coefficients 
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4. RESULTS 



The values for A, B, c\ \ and after tadpole 
improvement are shown in Fig. 1. The values for 
these parameters before tadpole improvement are 
typically very large. In the tadpole improvement, 
the perturbative expansion for the mean field pa- 

(2) (2) 

rameter «o = 1 + o^Wq , where u = —1.047, 
is applied. The renormalized coupling ay(q*) 
of Lepage and Mackenzie [8] is used to deduce 
values for a s ; fitting simulation measurements 
of the mean plaquette to its perturbative expan- 
sion specifies ay(q*) [5]. The computation of the 
scales q* for the above parameters is nearly com- 
pleted. 

It is possible to check the reliability of the per- 
turbative determination of the energy shift and 
mass renormalization. Recent NRQCD simula- 
tions [5-7] have produced values for the rest en- 
ergy aE^n(T) of the T. In addition, energies at 
small but non-zero momentum values have been 
obtained and fit to the dispersion relation 



#x(p) = #nr(T) 



2M- 
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yielding a nonperturbative estimate of the 

T mass. This mass can be compared with that 
expected from perturbation theory, namely, 



aM^ a) = 2(aM r - aE ) + aE NR (T), (2) 

where aEo = av(q*)A. The results of this com- 
parison are shown in Table 1. The masses agree 
down to the 1 — 2% level. 



5. CONCLUSION 

The leading radiative corrections to the heavy- 
quark mass renormalization, energy shift, and 
two important kinetic coupling coefficients were 
presented. The results underscore the reliability 
of tadpole-improved perturbation theory in the 
improvement of nonrelativistic lattice QCD. 
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